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ABSTRACT 

Context. Source I in the Orion-KL nebula is believed to be the nearest example of a massive star still in the main accretion phase. It is thus one 
of the best cases for studying the properties of massive protostars to constrain high-mass star formation theories. Near-infrared radiation from 
source I escapes through the cavity opened by the OMC1 outflow and is scattered by dust towards our line of sight. 

Aims. The reflected spectrum offers a unique possibility of observing the emission from the innermost regions of the system and probing the 
nature of source I and its immediate surroundings. 

Methods. We obtained moderately high spectral-resolution (A/AA ~ 9000) observations of the near infrared diffuse emission in several locations 
around source I/Orion-KL. We observed a widespread rich absorption line spectrum that we compare with cool stellar photospheres and 
protostellar accretion disk models. 

Results. The spectrum is broadly similar to strongly veiled, cool, low-gravity stellar photospheres in the range T e g ~ 3500 - 4500 K, luminosity 
class I-III. An exact match explaining all features has not been found, and a plausible explanation is that a range of different temperatures 
contribute to the observed absorption spectrum. The ID velocity dispersions implied by the absorption spectra, <r ~ 30km s _I , can be explained 
by the emission from a disk around a massive, m, ~ 10M Q , protostar that is accreting at a high rate, rri, ~ 3 x 10~ 3 M yr _1 . 
Conclusions. Our observations suggest that the near-infrared reflection spectrum observed in the Orion-KL region is produced close to source I 
and scattered to our line of sight in the OMC1 outflow cavity. The spectrum allows us to exclude source I being a very large, massive protostar 
rotating at breakup speed. We suggest that the absorption spectrum is produced in a disk surrounding a ~10 M Q protostar, accreting from its 
disk at a high rate of a few xl0~ 3 M Q /yr). 



1. Introduction 

Our understanding of massive star formation is still very lim- 
ited compared to that of solar-mass stars, for which a generally 
accepted theoretical framework exists ( |Shu et aLl [T987 ). The 
roles of dynamical interactions and coalescence versus a more 
standard core and disk accretion framework are still being de- 
bated (Beuther et al. 2007; Zinne cker & Yorke|2007[ >. Progress 
has been limited in part by the greater observational difficul- 
ties of studying massive protostars that are typically more dis- 
tant, embedded, crowded, and quickly evolving than their lower 
mass cousins. Orion is the closest star-forming region where 
young, massive (M> 8 Mq) stars are present, and as such is 
an ideal test case to study observationally and to understand 
theoretically. 

In particular, the Kleinman-Low (KL) nebula and the 
Becklin-Neugebauer (BN) object have received a lot of at- 



Neugebauer 1967 Kleinmann & Low 1967] >. Rieke et al. 



( 1973] > resolved several point sources in the KL nebula, and 



it was soon recognized that at least some of these objects were 
young massive stars. Subsequent higher angular resolution, in- 
frared observations showed that only a limited number of these 
sources were self-luminous objects. Infrared polarimetric ob- 
servations of the region have shown that one source is respon- 
sible for the illumination of most of the nebula, with the only 
exception being a small area around BN ( |Minchin et al.|199T I. 
Radio continuum observations have identified three sources 
that are likely young (possibly still forming) massive stars in 
the region: source n (also associated with an infrared source), 
BN, and source I, located near the Orion hot molecular core 
and near the infrared source IRc2 (M enten & Reid|19 95). The 
mid-IR spectrum of source n and its relatively low extinction 
imply that it is not particularly lu minous (|Gezari et al.|19 98). 
Similarly, the polarization data of Minchin et al. (1991) and 



tention since their discovery as infrared sources ( |Becklin &| [Werner et al.| (19831 indicate BN is not responsible for most 
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Based on observations collected at the European Southern 



Observatory, Chile. Program 076.C-0660. 



of the luminosity of the KL nebula. BN is known to be a run- 
away star (Plambeck et al. 1995), moving through the ONC 
at about 30 km/s. In one scenario, it was ejected from near the 
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Trapezium region by dynamical interaction with the 9 l C binary 
(Tan 2008 2004), and thus would have made a close passage 
to source I, perhaps enhancing its accretion rate by tidal inter- 
action. In an alternative scenario, BN was ejected by dynam- 
ical interaction with source I itself ( Bally & Zinnecker 2005 ; 
Go mez et al.|2 008 ). Source n may also have been involved in 



this interaction. This scenario requires source I to be composed 
of a massive, compact binary, with total mass more than twice 
the mass of BN, i.e. > 20 M Q . 



In either case, source I is thus expected to be the main illu- 
minating source of the infrared nebula. It is associated with an 
outflow detected in SiO and compact radio emission believed 
to originate in either a small ionized disk (Go ddi et al.| [2009; 
Plambec k et al.|19 90 2009 ) or in the ionized base of an out- 
flow cavity ( |Tan & McKee 2003). The source is deeply em- 
bedded and is not revealed at near infrared and shorter wave- 
lengths. Nevertheless, the surrounding reflection nebula offers 
a unique opportunity to observe the reflected emission from 
the heart of the system. [Morino et al. ( 1998 1 were the first to 
obtain a low spectral resolution, near infrared (2.05-2.35 //m) 
spectrum of the reflected light from source I. They revealed 
that the reflected spectrum shows the absorption lines typical 
of a cool photosphere, including metal lines and the CO(2-0) 
overtone absorption starting at 2.29yum. The inferred equiva- 
lent photospheric temperature (T+ ~ 4500 K) and derived lu- 
minosity (L* ~ 10 5 L ) could be consistent with either a very 
large, cool, massive protostar or the "photosphere" of a very 
active accretion disk. Nakan oet al.| ( [20 00) discus the first pos- 
sibility, using one-zone protostellar models and concluding that 
the required size of > 300/? o could not be achieved by rea- 
sonable accretion rates, even up to 10~ 2 M yr _1 . They sug- 
gest that the large size could be attained if the protostar was 
rotating at near the breakup speed, which would then imply 
very broadened line profiles. Subsequent work involving multi- 
zone protostellar models by Yorke & Bodenheimer (2008 ) and 
|Hosokawa & Omukai ( 2009 ) has shown that the sizes can be 
about a factor of several larger than those predicted by the one- 
zone models during the (relatively brief) luminosity wave ex- 
pansion phase (Stahler et al.|1986| of the evolution. jHosokawa 



& Omukai ( 2009 1 conclude that the low temperature implied 



by the infrared observations of the reflected spectrum from 
source I in Orion-KL could be achieved if the accretion rate 
was > 4 x 10~ 3 M yr~', averaged over the entire formation 
time of the protostar. 



However, such average accretion rates are at least an order 
of magnitude greater than those expected if massive stars form 
from massive cores in near equilibrium with the pressures in 



the Orion Nebula region (McKee & Tan 2003 1. Since a massive 



protostar will be surrounded by a relatively cool (compared to 
the protostar), optically thick accretion disk, we think that this 
is a more likely source of the observed near IR reflection neb- 
ula. The various models can be tested by resolving the pro- 
files of the absorption lines. With this in mind we carried out 
intermediate-resolution (R~ 8900), near-infrared observations 
of the reflected spectrum from source I in Orion-KL. 



2. Observations and data reduction 

Near-infrared medium resolution spectra of the scattered light 
around the Orion-KL region were obtained using the ISAAC 
near-infrared camera and spectrograph at the ESO-VLT UT1 
telescope. The observations were carried out in service mode 
on several occasions during October-November 2005. We used 
the 0.3 arcsec slit, and the medium-resolution grism that of- 
fered a resolution of ~8900 across the wavelength range 2.23- 
2.33 //m. The observing sequence was composed of sets of on- 
source and nearby (10 arcmin away) blank sky spectra arranged 
as standard ABBA cycles, to allow for efficient and accurate 
sky subtraction. For each of the two slit positions (see below 
and Fig. [TJ, we integrated a total of 40 min on source. Standard 
calibrations (flats and lamps) and telluric standards spectra 
were obtained for each observation as part of the ISAAC cali- 
bration plan. Wavelength calibration was performed using the 
lamp observations and refined using the bright OH sky lines. 

The spectra were reduced using standard procedures in 
IRAF. After checking that the data sets obtained in different 
nights were not significantly different, we combined all the ex- 
posures to obtain a final spectrum for each of the two slit posi- 
tions shown in the left panel of Fig. [T] 



3. Results 

Spectra were extracted for four regions corresponding to the 
brightest parts of the reflection nebulosity within the observed 
slits (see Fig. [T]for a detailed chart). Region "a" corresponds 
to the region observed by Morino et al. ( 1998} ; regions "b" and 
"c" correspond to two bright nebulosity patches to the west and 
to the south of the source I position, respectively; and region 
"d" is very close to the BN object. 

All spectra show emission lines of H2(2,1)S(1), 
H2(3,2)S(2) (blended with an unidentified infrared emis- 
sion feature), and [Felll]. These emission lines are most 



likely associated with the Orion HII region (Marconi et al. 
19981 |Walmsley et~al] [2000]) and/or the OMC1 molecular 



hydrogen outflow (e.g. 



Che rnin & Wrighf|[T996l |Kaifu et al. 



2000] |Nissen et al.| |2007| l. Some contamination from the 
molecular hydrogen emission can also be seen superposed on 
the CO overtone absorption. For the analysis of the reflected 
absorption spectrum, we have subtracted the nebular spectrum. 
We estimated the average nebular spectrum by averaging 
regions of the long slits away from the reflected emission and 
then subtracted it from each of the absorption spectra so as to 
remove the H2(3,2)S(2) line completely. 

The spectrum of region "d" is very different from the other 
three because it hardly shows any absorption feature, and the 
CO lines are detected in emission. This is consistent with the 
infrared spectra of the BN object (Scoville et al. 1979| >, and 
confirms that the nebular emission close to the BN object is 
dominated by the reflected light from this source rather than 
source I within the central KL region, as also inferred from 



the polarization maps of this region ( |Minchin et al. 19911. 
Therefore, we do not consider this spectrum further in the anal- 
ysis. 
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a (J2000) a (J2000) 

Fig. 1. Left panel: 2MASS Ks-band image of the Orion Trapezium and the Kleinmann-Low nebula region, 9 l C and the Becklin- 
Neugebauer object are marked for reference. The positions of the two on-source slits are shown. Right panel: the region around 
Orion-KL source I (marked with a cross and labeled) at 2.28/an is reproduced from an ESO-Archive SOFI-NTT observation 
(originally acquired for the ESO program 64.1-0493). The small slitlets marked from a to d show the regions we used to extract 
the spectra shown in Fig. [2] 



The polarization images of regions "a", "b", and "c" sug- 
gest that the emission here is dominated by scattered light 
from the embedded source or sources within the central KL re- 
gion (Mi nchin et al.|199l) . The three spectra (Fig. [2] top pan- 
els) show a common absorption line spectrum. In Fig. [2] we 
have marked the identification of the most prominent absorp- 
tion lines. In particular, along with the CO(2,0) and CO(3,l) 
ro-vibrational bands, we clearly detect the Cal triplet at ~ 
2.265 /an, the MgLl2.281 /an, and the FeLl2.239 /an lines. 

As mentioned above, our spectra from regions "a", "b", and 
"c" show a virtually identical absorption line pattern. This re- 



sult confirms and strengthens the conclusion by |Morino et al. 



( |1998| l that the absorption spectrum does not originate in the 
diffuse interstellar medium along the line of sight, but that 
we are instead observing the absorption spectrum of a single 
source that is scattered towards us by dust along different lines 
of sight towards the KL region. 

3.1. Comparison with cool stellar photospheres 



Mori no et al.| ( |1998| l points out that the reflected absorption 
spectrum they observed from the region roughly overlapping 
with our spectrum "a" has an absorption line pattern similar to 
that of a cool stellar photosphere. Their analysis suggests that 
the absorption spectrum, especially the atomic lines, is consis- 
tent with that of the giant star HD12014 veiled by an excess 
featureless continuum emission, which contributes about half 
of the continuum emission at ~2 fim. HD12014 has a KOIb pho- 
tosphere, corresponding to an effective temperature of ~4500 K 
and low gravity. 

Our higher resolution spectra allow a more detailed com- 
parison with stellar photosphere spectral libraries. In Fig. [3] 



we show the comparison between the spectra in the library of 
Kleinmann & Hall ( 1986| > and our average reflected spectrum, 
obtained by combining the data from regions "a", "b", and "c". 
Our spectrum was smoothed to the same resolution as the spec- 
tra in the library. The photospheric spectra were veiled in such a 
way to approximately match the depth of the central line of the 
Ca I triplet in the Orion spectrum. To match the Ca I absorption 
the amount of veiling required by the late type spectra is such 
that the stellar photosphere contributes ~ 20% of the computed 
spectrum, while for the earlier spectral types the photosphere 
dominates the spectrum. 

The ratio of the atomic absorption lines is best matched 
with veiled stellar photospheres of spectral types M2-K0, cor- 
responding to effective temperatures in the range 3500-4500 K. 
Later types have much more pronouced Cal absorption than 
the Mgl, while the opposite is true for earlier spectral types. In 
all cases, as noted by Morino et al. ( 1998| l, the CO absorption 
bands are much more pronounced in the Orion spectrum than 
in the stellar photospheres. In this respect, the lowest gravity, 
later type spectra in the library are closer to our observed spec- 
trum. This suggests that the spectrum may be consistent with 
an even lower gravity "photosphere", such as from an accretion 
disk. 



3.2. Velocity dispersion of the absorption spectra 

To estimate the broadening of the absorption line spectra, we 
compared the CO absorption in our ISAAC spectra with the CO 



absorption in the high-resolution spectra of Wallace & Hinkle 
(1996). The spectra from the library were convolved with a 
Gaussian of variable width to mimic line broadening (in ad- 
dition to that intrinsic to the stellar library), shifted in velocity, 
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Fig. 2. From top to bottom: final spectra of the four regions a to d shown in Fig. [T] Each spectrum has been normalized to the 
continuum level. Some of the most prominent emission or absorption features are labeled. The features labeled "Telluric" are 
negative residuals from the subtraction of the OH airglow spectrum. The spectra also show the prominent emission features of 
H2 and [Felll] which originate from diffuse gas in the region and (most likely) not associated with the reflection nebulosity. 



veiled and then smoothed at the same resolution and resampled 
onto the same wavelength grid of our ISAAC spectra for proper 
comparison. These parameters were allowed to vary and ax 2 
cube was constructed. In Fig.|4]we show, for the fit to the spec- 
trum of a Ori, the variation of the^ 2 as a function of the ID ve- 
locity dispersion (cr = FWHM/ 1.665) of the Gaussian smooth- 
ing function for each of the three regions "a", "b", and "c". The 
results obtained using this star are consistent with those of the 
other late type, low-gravity stars in the library (such as a Her or 
RX Boo). The results show that to match our spectra, the pho- 
tospheric templates need to be smoothed with Gaussian profiles 
with cr < 30 km/s for region "c", while the spectra of regions 
"a" and "b" require a much smaller smoothing, cr < 15 km/s. 
The true, total velocity dispersions will be greater than these 
limits owing to the contribution of the intrinsic dispersion in 
the stellar atmospheres used in the spectral library. 

Another estimate of the line broadening can be done di- 
rectly with a simple modeling of the CO lines. The line param- 
eters were taken from the HITRAN database (Rothman et al. 
2009), and line intensities were computed following Helmich 
(1996) using a simple LTE approximation. As in the proce- 
dure with the template spectra, we convolved the model spectra 
with a Gaussian profile of a given full width at half maximum, 
and the synthetic spectrum was then scaled to match the ob- 
served absorption spectrum, smoothed at the appropriate spec- 
tral resolution, and resampled on the same wavelength grid as 
our ISAAC spectra. This simplified analysis does not allow us 
to accurately reproduce all the CO line intensities, but the line 



broadening required to match the observed spectra is consistent 
with the results derived using the spectral templates. In partic- 
ular region "a" requires cr < 20 km/s, while "b" and "c" are 
consistent with cr ~ 27 + 3 km/s. 

We repeated the analysis using a double-peaked line profile 
as expected from a disk (e.g. Kenyon et al. 1988 ). Although we 
do not observe double-peaked lines, our observations would 
still be consistent with line profiles where the double peaks are 
separated by < 45 km/s. 

The estimates of the linewidths are also consistent with re- 
cent observations of mid infrared CO absorption lines toward a 
region close to our region "b" (Beuther et a l.|2010 their source 
IRC3). Their spectrum also shows emission components, most 
likely from outflowing gas, which further complicates the anal- 
ysis of the line profile. 



4. Discussion 

Our higher resolution spectra (R~8900 vs. ~1000) allow us 
to test the different scenarios p ut forward by | Morino et al.| 
( fT998) l, |Nakano et al.| ( |2000l ), and |Hosokawa & Omukai|p009| l 
on the nature of the source producing the absorption spectrum. 
The latest authors have attempted to explain the observed, cool 
spectrum as being produced by a very large (r, > 200 R Q ), 
hence cool (T < 5500 K), protostar, which requires very 
high accretion rates averaged over the whole formation time 
(m* > 4 x 10~ 3 M yr~'), and the protostar to be at a particular 
stage in its evolution (m, ^ 2OM ), undergoing swelling due to 
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Fig. 3. Comparison between our average reflected spectrum (solid line) and the stellar spectra from 
dotted lines). In the left panels we show enlargement around the lines of Fe I, Ca I, and Mg I (the wavelength of each feature is 
marked with a vertical green dotted line), and on the right side we show the full wavelength range of our spectrum, and the name 
and spectral type of the stars are marked in the right side panels. 



the "luminosity wave" (Stah leretal.|19 86) immediately before 
the onset of Kelvin-Helmholz contraction. 

For a protostar forming from a near-equilibrium gas core, 
such high average accretion rates would require very high am- 
bient pressures, equivalent to the self-gravitating w eight of a 
gas clump with E > 44 g cm" 2 (McKee & Tan 2003 1. However, 
as described by these authors, the Orion Nebula Cluster has a 
much lower mean mass surface density of 0.24 g cm" 2 . While 
moderate enhancements in £, pressure, and thus m* are ex- 
pected in the central regions of the ONC, it appears to be diffi- 
cult to create the conditions that would allow such high accre- 
tion rates to produce the observed NIR scattered light directly 
from the protostellar surface. The instantaneous accretion rate 
from a turbulent core is expected to show significant fluctua- 
tions about the mean value. Fluctuations may also be induced 



by tidal interactions with passing stars in the cluster. This sce- 
nario is discussed below in the context of source I and the BN 
object. However, these fluctuations are not expected to signifi- 
cantly change the size of the star, which is set by the accretion 
rate averaged over the last growth time (m»/m») of the proto- 
star. 

The large protostar scenario makes predictions for the ve- 
locity dispersion of the CO banheads as a function of the rota- 
tional velocity of the protostellar surface relative to its breakup 
velocity. For example a 20 M Q , 200 R protostar rotating at 
50% of the Keplerian velocity at its surface has an equatorial 
rotation speed of 69 km s _1 . Viewed at an inclination angle of 
45°, this would produces an observed line profile with a veloc- 
ity dispersion of a tat =s 50km s" 1 (e.g. Gray 19921. 
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Fig. 4 . Results of the^- 2 fits to the CO absorption in Orion with the high-resolution infrared spectrum of a Ori (Wallace & Hinkle 
1996 1. In the left side panel we show the Ax 2 curves as a function of the AV used to smooth the library spectrum. On the right 



we show the comparison between the observed spectrum in each of the three regions (full lines) with the best-fitting photosphere 
(dashed lines, corresponding to the A^ 2 minima in the left-hand side plot). 



While the large protostar scenario is still marginally con- 
sistent with the observed velocity dispersion, we favor the al- 
ternative scenario of the NIR spectrum showing CO bandhead 
absorption being produced in an active accretion disk. There 
are several arguments in support of this: (1) accretion disks are 
expected to be present around protostars; (2) the large protostar 
scenario requires very high accretion rates and a very special 
protostellar evolutionary phase; (3) the detection of prominent 
low excitation CO absorption near 4.6 /im in the reflected spec- 
trum from region "b" (IRC3) by Beuther et al. (2010) suggests 
there is a range of "photospheric" temperatures. This is ex- 
pected in an active disk, but difficult to reconcile with a single 
stellar photospheric temperature — in fact, these low excita- 
tion lines are not prominent in cool stellar photospheres ( Sloan 
|et al.|2003] e.g.); (4) there is evidence that source I is in an ac- 
tive accretion phase, such as the intense outflow activity from 
the region (e.g. the OMC1 outflow) and the claims that the cm 
radio emission and SiO maser motions around source I indicate 
the presence of either a compact circumstellar disk (Reid et al 



1996 1 and the morphology of the NIR nebula. We suspect there 
may still be the possibility of misinterpreting of the outflow and 
disk orientations from the maser data since the maser emission 
may not provide a uniform sampling of the gas in these struc- 
tures. Alternatively, if there has been a strong, recent dynamical 
interaction of source I with BN ( B ally Ik Zinnecker|2005 1 Tan 
2004), then the orientation of the disk may have been changed. 
However, it is still difficult to understand how direct illumi- 
nation of the NW and SE regions of the KL nebula could be 
achieved by source I in this scenario, so for these reasons we 
favor the fiducial model as we have described it. 

4.1. Comparison with the turbulent core protostellar 
model 



2007 1 Matthe ws et al.|20"T0) > or the ionized base of an outflow 
cavity (Tan & McKee 2003). 

The accretion disk scenario was also briefly discussed by 



McKee & Tan ( 2003 \ have estimated the range of protostellar 
models consistent with producing the observed luminosity of 
the Orion Hot Core, i.e. the Kleinmann-Low Nebula, of ~ 1 - 
5 x 10 4 L Q . Here we consider a comparable range of models, 
with parameters listed in Table [T] They involve a 60 M gas 
core forming a star with 50% efficiency and at various stages 
of collapse, i.e. with m, in the range of 8 to 20 M Q . If the core 
Nakano et aL] ( |2000| l. To further explore this possibility, in the was in pressure equilibrium with a I = 1 g cirT 2 surrounding 

clump, then the accretion rates would be a few xl0~ 4 M yr _1 . 

It is possible that the Orion KL protostar has suffered a tidal 
perturbation in the past ~ 1000 yr due to the close passage of 
the runaway BN star ( Tan|20 04). This is expected to have en- 
hanced the accretion rate through the inner disk and thus also 
the protostellar mass outflow rate, perhaps causing the appar- 
ently "explosive" morpholgy of the outflow (Allen & Burton 



next section we compare our spectrum with a more detailed 
model of a massive protostar-disk system as expected in the 
turbulent core model for high-mass star formation (Mc Kee &| 
|Tan|2003D . 

There are claims that the outflow and disk from source I are 
oriented orthogonally to the orientations assumed in our model, 



i.e. with the outflow aligned with a NE-SW axis (e.g. Goddi 
|et al. 2009, Reid et al. 2007). Such an orientation is hard to rec- 
oncile with the observed orientation of the large-scale outflow 
from the KL region ( |Allen & Burton|1993[|Chernin & Wright 



1993 ). Thus we also consider some models with accretion rates 
enhanced by a factor of 10, especially for lower mass protostars 
that will still satisfy the observed total luminosity constraint. 
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For the 15 M G protostar accreting at m* 



3.28 x 



10 M Q yr , the star has a mean spherical radius of 7. 4-2R Q and 
internal luminosity of L, = 2.01 x 10 4 L Q with a photospheric 



temperature of 25,400 K ( |McKee & Tan|2003[ ). The accretion 
boundary layer has a luminosity of Lbl — 1-03 x 10 4 L Q (with 
a photospheric temperature of 32, 000 K). The accretion disk 
is assumed to have the same luminosity as that of the bound- 
ary layer. These estimates of accretion luminosity should be re- 
garded as upper limits since they do not include any reduction 
from the powering of protostellar outflows, which may have 
mechanical luminosities of about 50% of the total accretion 
power. 

We assume the protostar has a geometrically thin, opti- 
cally thick accretion disk. We consider the case of a purely 
active disk, i.e. with negligible heating of the disk by the pro- 
tostar, and the case of a protostellar-heated disk, following the 
treatment of |Matzner & Levin] ( |2005| > and |Kratter & Matzner| 
(2006). The reprocessing factor, f{ n - is defined as the ratio be- 
tween the incident flux F ln normal to the disk surface, and the 
spherical stellar flux at that radius: F m = f^L/iAur 2 ). An an- 
alytic approximation for f m is f m 0.1 , where e is the 
instantaeous star formation efficiency from the core, which we 
have taken to be 50% ( |Matzner & McKee||2000| l. In both the 
non-irradiated and irradiated cases, emission from the accre- 
tion disk dominates the emergent NIR flux, as shown in Fig. [5] 

The radial profiles of disk photospheric temperature and 
Keplerian velocity are shown in Fig. [6] We evaluated the 
contribution to the total CO line profile from each emitting 
part of the disk, weighting by the disk area, flux at 2.3 yum, 
and a temperature-dependent "CO bandhead weighting factor", 
fco(T), which is derived from the observed equivalent widths 
of giant and supergiant standard stars (see Fig. FT}. This factor 
is uncertain at T < 2000 K, although we have checked various 
functional forms that decline with temperature for this temper- 
ature range and found that our results are not particularly sen- 
sitive to such variations. We assume a typical inclination angle 
of 45° of dust particles to the rotation axis of the accretion disk, 
since the RA. of the axis of the bipolar flow (and the observed 
near IR reflection nebulosity) is approximately from the NW 
to the SE (e.g. |Chernin & Wright|[T996| . The application of 
this outflow geometry to source I is consistent with the inter- 
pretation of source I's radio emission as the ionized base of an 
outflow cavity (Tan & McKee 2003), but not with the interpre- 
tation that the disk axis is oriented from the NW to SE (Reid 
et al.|2007[ [Matthews et al.pjTD) . 

The ID velocity dispersions, <x, of the resulting line profiles 
are listed in Table Q] The irradiated disks have lower values 
caused by the disk being warmer at a given radius and the CO 
emitting region being located at larger radii. If an inclination 
angle of 30° is adopted, for example, for the 15 M model with 
m„ = 3.28 x 10~ 4 M Q yr _1 , then the velocity dispsersion drops 
to <r = 49.1 km s _1 and cr ilT = 40.1 km s _1 . 

The models with enhanced accretion rates cause the disk 
to be warmer at a given radius, thus pushing the CO bandhead 
region of the disk to larger radii and lower Keplerian velocities. 
This produces lower values of cr. 

We conclude that the fiducial turbulent core protostellar 
models (i.e. those without enhanced accretion rates) for the 



Orion KL protostar predict CO bandhead ID observed veloc- 
ity dispersions that are in the range 53 to 75 km s _1 . The ve- 
locity dispersion is not very sensitive to m* because the tem- 
perature range that contributes to CO bandhead absorption, 
i.e. T ~ 2000 - 4000 K occurs at a fairly constant range of 
Keplerian velocities in the accretion disk, nearly independent 
of ra,. Disk irradiation raises the temperature of the disk at 
any given radius, thus shifting the zone that contributes to CO 
bandhead equivalent width to larger radii and thus yielding nar- 
rower line profiles. The velocity dispersion of the line profile is 
somewhat sensitive to the adopted inclination angle between 
the scattering surfaces in the outflow and the disk rotation axis. 
The fiducial protostellar models predict ID CO bandhead ve- 
locity dispersions that are about a factor of two to three greater 
than observed. A better match is achieved if there is a relatively 
stronger contribution to CO bandhead EW from lower velocity 
dispersion gas, which can be reproduced in models with accre- 
tion rates enhanced by a factor of 10. Higher velocity resolution 
observations of the scattered light from Orion KL are needed 
to better establish the significance of the observed velocity dis- 
persion in order to further test theoretical models, including the 
search for double-peaked line profiles. 



5. Summary and conclusions 

In this paper we have presented medium-resolution (R~ 8900), 
near-infrared spectra of the reflection nebulosity in Orion-KL. 
Our main results are as follows. 



- The absorption line spectra in various regions of the re- 
flection nebula in Orion-KL are consistent with each other, 
except for a small difference in the line broadening. This 
suggests a common origin for the absorption line close to 
the region emitting the continuum emission. We thus con- 



firm the earlier suggestion by Morino et al. ( 1998 1 that the 
absorption spectrum is generated close to the position of 
source I, either from the protostellar photosphere or the sur- 
rounding accretion disk. This also suggests that the outflow 
cavity from the main illuminating source in the KL nebula, 
i.e. source I, is aligned from the SE to the NW, which is 
difficult to explain in models that invoke a disk plane along 
this axis ( |Reid et al.|2007)|Matthews et al.|20T0) . 
- The observed absorption spectra are broadly consistent 
with low-gravity, late-type, heavily veiled stellar photo- 
spheres with spectral types in the range M2-K0 and lumi- 
nosity class I-III, corresponding to effective temperatures in 
the range 3500-4500 K. Nevertheless, it was not possible to 
find a good match for all the observed absorption features 
with a single photospheric template. In particular, the ob- 
served CO absorption in the Orion spectra are much deeper 
than in the stellar photospheres that match the atomic lines. 
Although it is possible that an even lower gravity photo- 
sphere than those used in the comparison may alleviate or 
solve this inconsistency, it is likely that this discrepancy 
could be solved if the Orion spectrum is produced by the 
combination of a wide range of photospheric temperatures, 
as in a very active accretion disk. 
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Table 1. Parameters of protostellar models 
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Fig. 5. Spectra of protostellar models (see text). Top panel: 
For the fiducial accretion rate model with m, = 10 M , the 
fluxes from the star and boundary layer are shown by the 
dashed and lower long-dashed lines, respectively. The lower 
solid line shows the flux from the accretion disk (integrating 
out to 5000r») assuming negligible protostellar heating, while 
the lower dotted line shows the case accounting for proto- 
stellar heating. The upper long-dashed, solid, and dotted lines 
show the equivalent quantities for the model with the accre- 
tion rate enhanced by a factor of 10. Bottom panel: The solid 
lines show for the fiducial m* = 10 M Q model the flux from 
the star -(-boundary layer, nonirradiated and irradiated accre- 
tion disks. The long-dashed, dotted, and dashed lines show the 
equivalent fluxes from the fiducial accretion rate models with 
m* = 8, 15,and20 M Q , respectively. In all cases, the flux at 
2.3 fim is dominated by disk emission. 



- The observed line profiles are consistent with being un- 
resolved at our spectral resolution in regions "a" and "b" 
and marginally resolved in region "c". These results im- 
ply cr ~ 30 km/s for region "c" and about a factor of two 
less for regions "a" and "b". We thus exclude that the spec- 
tra could be originated from the photosphere of a massive 



Fig. 6. Radial profiles of protostellar accretion disk models. 
Top panel: Disk photospheric temperature for non-irradiated 
disk (lower solid line) and irradiated disk (upper solid line) 
for m, = 10 M with fiducial accretion rate. The equivalent 
models for m, = 10 M Q with accretion rate enhanced by a fac- 
tor of 10 are shown by the dotted lines. The equivalent models 
for m, = 20 M Q with fiducial accretion rate are shown by the 
dashed lines. The long-dashed horizontal lines show the ap- 
proximate range of temperatures where we consider there are 
contributions to the CO bandhead equivalent width. Bottom 
panel: Keplerian velocities for the m* = 10,20 M models, 
shown by the solid and dashed lines, respectively. 

(~ 20M Q ), large (~ 200/? o ) protostar rotating close to the 
breakup speed. The linewidths are consistent with the pre- 
dictions of a model of an active accretion disk around a 
massive protostar (M~ 10 M ). Models with accretion rates 
enhanced above the fiducial values are favored in order to 
reproduce the relatively narrow line widths. 

We conclude that the most likely source of the near infrared 
absorption spectrum is the accretion disk surrounding the mas- 
sive protostar source I, with the system in a state of relatively 
high accretion rate. 
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Fig. 7. Weighting factor of CO(2-0) bandhead equivalent width 
over the wavelength range 2.292 - 2.322 //m derived from gi- 
ant/supergiant standard spectra (open squares) and Mira spectra 
scaled by a factor of 0.7 (solid squares). 
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